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O^ ; ABSTRACT 

We consider time-periodic nonlinear localized excitations (NLEs) on one- 
dimensional translationally invariant Hamiltonian lattices with arbitrary fi- 
Q I nite interaction range and arbitrary finite number of degrees of freedom per 

Y] I unit cell. We analyse a mapping of the Fourier coefficients of the NLE so- 

lution. NLEs correspond to homoclinic points in the phase space of this 
^ ■ map. Using dimensionality properties of separatrix manifolds of the map- 

ping we show the persistence of NLE solutions under perturbations of the 
system, provided NLEs exist for the given system. For a class of noninte- 
grable Fermi- Past a- Ulam chains we rigorously prove the existence of NLE 
solutions. 
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1 Introduction 

The existence, stability and properties of nonlinear localized excitations (NLEs) 
in Hamiltonian lattices with discrete translational symmetry have been sub- 
jects of growing research interest (see e.g. |jl| , , ||^ , and references therein) . 
NLEs can be viewed as generalized discrete analogs to the breather solution 
in a sine-Gordon equation 0. They are characterized by a localized vibra- 
tional state of the lattice. There are two basic reasons for the generic NLE 
existence on Hamiltonian lattices: i) the lattice forces acting on a given par- 
ticle are nonlinear (thus one can tune oscillation frequencies by varying the 
energy) and ii) the discrete translational symmetry of the lattice (in contrast 
to the continuous translational symmetry of Hamiltonian field equations) 
provides a finite upper phonon band edge of the spectrum of extended small 
amplitude oscillations of the lattice around its groundstate [§],0|. 

There are a few known rigorous NLE existence proofs. First NLEs are 
exact solutions of the integrable Ablowitz-Ladik lattice 0]. In fact they form 
a three-parameter family of solutions. Secondly NLEs are exact solutions 
for the Fermi-Pasta-Ulam chain with box-like interaction potential |^. In 
this case the NLEs are of compact support. Finally, and most importantly, 
MacKay and Aubry have derived an existence proof for NLEs in an array of 
weakly coupled anharmonic oscillators 0. Remarkably this existence proof 
works independent of the lattice dimension. 

In this contribution we will first deal with the existence of NLEs in nonin- 
tegrable generic one-dimensional Hamiltonian lattices. In the second part we 
will investigate a class of Fermi-Pasta-Ulam chains and give rigorous proofs 
for the NLE existence. 

Let us briefly outline the main steps in the approach presented below. 
We will assume the existence of a time-periodic NLE on a one-dimensional 
lattice. We represent the lattice displacements at each lattice site in the 
NLE ansatz in a Fourier series with respect to time. When we insert the 
NLE ansatz into the lattice equations of motion, we obtain a set of coupled 
algebraic equations for the Fourier components of the NLE ansatz, which 
form an infinite-dimensional map. The NLE solution has to correspond to a 
common point of two separatrix manifolds in the phase space of the map, or 
a homoclinic point. Analyzing the map in the (linearizable) tails of the NLE, 
we can derive the dimension of the separatrix manifolds. Consequently we 
show that if a homoclinic point exists for a given system, then generically the 



hoinoclinic point will survive under perturbations of the system. We then 
consider a subclass of Hamiltonian chains and rigorously prove the existence 
of two different (with respect to symmetry) NLE solutions. For this particular 
example we show the emergence of horseshoe patterns - a consequence of the 
existence of homoclinic points. 

2 Stability of NLE solutions under Hamilto- 
nian perturbations 

We consider a classical one- dimensional Hamiltonian lattice of interacting 
particles (perhaps feeling an external field periodic with the lattice) with 
lattice cite a = 1. The displacements of the particles from their groundstate 
(equilibrium) positions are given by a ra-dimensional vector Xi, where n is 
the number of components per unit cell {n < no, uq finite) and the integer I 
marks the number of the unit cell. The range of the interaction r is consid- 
ered to be finite: r < tq, tq finite. Here no and tq are positive integers. The 
potential energy of the system is required to have anharmonic terms in the 
displacements if expanded in a Taylor series around the groundstate (mini- 
mum of potential energy) of the system. Furthermore the potential energy 
should become a positive definite quadratic form in the limit of infinitely 
small displacements. The Hamilton function H is given by the sum of the 
kinetic energy of all particles and the potential energy. 

As it was shown in 0,0' ^^^ only possible exact NLE solution on an 
arbitrary lattice (by that we mean no additional symmetries are present) has 
to have the form: 

Xi{t) = Xi{t + 27r/ui) , Xi^ioo ^ . (1) 

Then one can avoid resonance conditions of multiples of the fundamental 
frequency ui (as they appear in the Fourier transformed functions in (|I|) 
with respect to time because of the nonlinearity of the system) with phonon 
frequencies of the linearized (around the groundstate) system. Since the 
motion of assumed existent NLEs requires the excitation of at least a second 
fundamental frequency in the ansatz (|l|) JTOf, we can exclude them from the 



consideration and search for stationary time-periodic NLEs as given in (|1|). 
Note that in the case of the Ablowitz-Ladik lattice additional symmetries 



are present (the lattice is integrable) and thus the above statement does not 
hold - moving NLEs are exact solutions in this nongeneric case. 

Because of the assumed time-periodicity of all displacements in (|I|) we 
obtain for the assumed NLE solution 

^1= E Ai,k)e'^''' ■ (2) 

A:=0,±l,±2,... 

Now we can insert this ansatz (|^) into the Newtonian equations of motion 

Xi = dH/dXi. The left and right hand expressions of the equation of motion 
are represented again as a general Fourier series. Equaling the prefactors at 
identical exponential terms to each other we finally obtain a coupled set of 
algebraic equations for the unknown Fourier coefficients A{i,fc)- Because of (|^) 
the Fourier coefficients have to satisfy the boundary condition A(^i^±oo,k) — ^ 0. 
In the following we will study properties of this algebraic set of equations. 

First we can note that because of the d = 1 dimensionality of the consid- 
ered lattice the coupled set of algebraic equations for the Fourier coefficients 
y4(/ fc) can be represented as a discrete map of a cIm = (2ronofcmax)-dimensional 
phase space, where the integer kmax represents the number of considered 
higher harmonics in (0) and has to tend to infinity. In other words, given 
the Fourier coefficients at 2ro neighboring lattice sites completely determines 
the Fourier coefficients to the left and right of the specified chain segment. 
Secondly because of the required asymptotic vanishing of the Fourier coeffi- 
cients for / — > ±oo we can linearize the dM-dimensional map in the tails of 
the assumed NLE solution with respect to the variables ^(/,fc)- The linearized 
map will decouple into a set of kmax independent ds = 2rono-dimensional lin- 
ear submaps p. In each of these submaps Fourier components with only 
one Fourier integer k will appear. Each linear submap is equivalent to the 
problem of finding solutions of the linearized (around the groundstate) lat- 
tice equations of motion using the ansatz Xi{t) = Aiexpiujt. Here for every 
submap one has to substitute u = kui. If oj equals to an arbitrary phonon 
frequency of the linearized lattice equations, then it follows that the corre- 
sponding linear submap is characterized by a matrix Mfc of dimension ds x ds 
with detMfe = 1 and AjAj+rf^ = 1 where Xi,i+ds are two of the 2ds eigenvalues 
of Mfc. Since these properties do not depend on the specific value of the con- 
sidered phonon frequency, it follows that they are independent on uj and thus 
independent on k and uji. Consequently we find that every linear submap is 
volume preserving and symplectic. 



Finally we note that for every considered linear submap (and thus also 
for the original (iM-dimensional map) the phase space point of zero Fourier 
components is a fixed point of the map. If kui equals a phonon frequency, 
the fixed point is of elliptic character. If however kui does not equal any 
phonon frequency, the fixed point has to be a saddle point (because it can 
not be an elliptic fixed point and the map is symplectic). That means that 
tqUq eigenvalues of the matrix M^ are real and of absolute value lower than 
one. Consequently the separatrix manifold of all points in the space of the 
submap attracted by the saddle point is of dimension dsm = roUQ. 

If we require that neither of the multiples of the fundamental frequency 
kuji resonates with the phonon band, the corresponding original ^M-dimensional 
map in the phase space of the Fourier coefficients has a separatrix manifold of 
dimension c^a/ /2. To get a NLE solution we have to find a point in the space 
of the original map which belongs simultaneously to the separatrix manifold 
S^ attracting the solution to zero for I -^ — oo and to the separatrix mani- 
fold 5+ attracting the solution to zero for / -^ +00 (cf e.g. ||ll|,[|l3)- This 



point has to be then a homoclinic point |T^. If a homoclinic point exists, 
there exist an infinite number of homoclinic points, which can be obtained by 
subsequent mapping of a given homoclinic point. As a result the horseshoe 
structure of intersections between the stable and unstable manifolds have to 



13 



emerge 

Because of the assumed nonresonance condition (see above) every mani- 
fold has dimension (iM/2. Let us choose one of the homoclinic points. Then 
we can define the tangent planes to each of the two manifolds in this point. 
There are two possibilities for the topology of these two planes. Either i) 
these two planes span the whole phase space of the map (of dimension c?m) 
or ii) the two planes span a space of lower dimension than dM- In case i) 
any small perturbation of the original system (consequently of the map, con- 
sequently of the two manifolds and consequently of the two tangent planes) 
will only shift the homoclinic point smoothly. Also in this case i) all other 
homoclinic points will have the same topology with respect to the tangent 
planes. In case ii) there exist perturbations of the system which will lead to 
a vanishing of the homoclinic point (and thus of all other homoclinic points 
too). Still there will exist a perturbation for case ii) such that the homoclinic 
point is smoothly shifted and simultaneously the two tangent maps will span 
the whole phase space of dimension dM- 

Consequently we can conclude, that if we have a NLE solution which cor- 



responds to a set of homoclinic points of type i) then any small perturbation 
of the system will smoothly transform the NLE solution into a new NLE 
solution. If we have a NLE solution which corresponds to a set of homoclinic 
points of type ii), we can always find a perturbation such that we transform 
the NLE solution into a new NLE solution which corresponds to a set of 
homoclinic points of type i) and is then stable under subsequent small per- 
turbations. From the above said it follows that NLE solutions, if they appear, 
are generically not isolated objects in the sense that small perturbations of 
the system either do not destroy them at all, or that a proper perturbation 
of the system will transform them into nonisolated solutions. 

We have operated with finite phase space dimensions dM- Of course (Im 
was not bounded from above, so that the limit (Im — >■ C)0 can be considered 
without altering the arguments. 

Let us summarize the results obtained so far. If we assume that for 
a given Hamiltonian chain a NLE solution exists, then it corresponds to 
an infinite number of homoclinic points in the map as introduced above. 
If neither of the multiples of the frequency of the NLE solution resonate 
with the phonon band (nonresonance condition), then the dimension of the 
stable and unstable manifolds is exactly one half of the map's phase space. 
From topological arguments it follows then, that either the NLE solution is 
structurally stable - i.e. it is smoothly transformed under any perturbation 
of the Hamiltonian of the system, or there exists at least one perturbation of 
the Hamiltonian such, that the NLE solution is smoothly transformed into 
a structural stable one. Of course a perturbation of the Hamiltonian has to 
preserve the general structure of the map, i.e. we are not allowed to consider 
e.g. a two-dimensional perturbation, which would lead to the impossibility 
of defining the map. 

A central point in the consideration so far has been the nonresonance 
condition. If this condition is violated for n values of kui, then the dimen- 
sion of the stable and unstable manifolds is lowered by {2n) to {dnj — 2n). 
If an NLE solution exists for such a case, then the tangent planes of the two 
manifolds in a given homoclinic point can not span the whole phase space of 
the map. Consequently there will be always perturbations of the Hamilto- 
nian such that the NLE solution will vanish. At this point it also becomes 
clear that in the analogous problem of a Hamiltonian field equation, where 
resonances can be never avoided [|12| , systems which allow for NLE solutions 



become isolated (see also jH^ ). 
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If a system is given with a NLE solution which fulfills the nonresonance 
condition, it could be possible that a given perturbation of the Hamiltonian 
would lead to a violation of the nonresonance condition. Consequently one 
has to check in a given case, whether the nonresonance condition survives 
under a perturbation. 

3 Proof of existence of NLE solutions for Fermi- 
Pasta-Ulam chains 

In the last part of this work we will prove rigorously the existence of NLE 
solutions for a class of Fermi-Pasta-Ulam chains governed by the following 
equations of motion: 

Xi = -{Xi - X,_i)"""' - (Xi - X,+i)'™"' (3) 

with m = 2,3,4 As it was shown in 101,11131 we can consider a time- 



space separation ansatz Xi{t) = AiG{t). Inserting the separation ansatz 
into (|^) we get a differential equation for G(t): G(t) = — fi;G^™~^(t) and a 
two-dimensional map for the amplitudes Af. 

nAi = (A - A,_i)'™"' + (A - A+,Y^-' . (4) 

Here k > is required in order to get a bound oscillatory solution for G{t). 
The phase space properties of map (^ are shown in Fig. 1 for the case m = 2 
and K = 1. Below we will refer to particular patterns observed in Fig.l. We 
consider cases when Ai/Ai^i < and introduce fi = \— l\^Ai with 

«:// = (// + //-i)'"^' + (// + fi+i)'""-' . (5) 

Equation (^) can be viewed as a two-dimensional map M of a vector // = 
if I, fi-i)'- fi+i = M/;. The task is then to show the existence of at least one 
value of K > such that (||) yields f\i\ > f\i/\ > for |Z| < |/'| and fi^±oQ —>■ 0. 
First we note that for any value of k, the map (|^) has a fixed point 
ffpi = (0,0). Adding weak harmonic nearest neighbour interactions to (|^) 
and considering the limit of vanishing harmonic interactions yields that the 
fixed point (0, 0) is a saddle point with eigenvalues Ai jpi = and X2,fpi = 
cxD. Consequently we find that there exist two one-dimensional separatrix 
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manifolds 5"+ and S- of the map (^. All points in the two-dimensional 
phase space of the map belonging to S± are attracted to the saddle point 
after an infinite number of iterations for / -^ ±00. Secondly it follows that for 
Kfp2 = 2^*"/^™"^ the point ffp2 = (/, /) is also a fixed point (see also Fig.l). 
Linearizing the map around ffp2 yields the elliptic character of ffp2. The two 
eigenvalues of the linearized map Aijp2 and X2jp2 obey the relations |Aijp2| = 
|A2,/p2| = 1 and are given by the expression Ai2,/p2 = h — l±i{l — (1 — k)^)^/^ 
with k = 2/(2m- 1). 

Let us mention another useful property of the map (|), which holds due 
to the discrete translational symmetry of the considered system @. If we 
generate a sequence of /2, /a, ... starting with a vector /i = (a, b) and iterating 
(^ towards growing integers /, then we can generate the same sequence of 
numbers by iterating (||) towards decreasing integers I starting with the vector 
fi = {b,a). 

We consider an initial vector fi = {a, h) with < a < 6 and iterate 
towards increasing /. If we require < //+i < a, the parameter n is confined 
to Ki <K< 4 With K^ = {{fi + fi-if^-' + /,''"-^)//z and 4 = kT + 
^22m-i _ i)y2m 2^ p^^ ^^^ ^£ ^-^^ obtained values of Z^+i = c we can consider 
the next step requiring < //+2 < c. By itself this requirement again yields 
a confinement of k to n^+i ^ '^ ^ '^t+i- ^^ ^^ choose c = then k~[j^^ = +cxo. 
If we choose c = a it follows 4+i — '^t ■ Because c is a monotonic function 
of K we can satisfy a > //+i > //_|_2 > only if kJ < n < 4 for a = b and 
nj < n < 4 for a < b. Repeating this procedure for every following lattice 
site we get a sequence of monotonically increasing lower bounds on n defined 
by the vanishing of fi+m = in the m-th step. In the limit m ^ oo we thus 
obtain exactly one value of k, which is the limit of the mentioned sequence of 
lower bounds. This limiting value is finite, because it is always smaller than 
4. It can not be equal to 4 either, since we would then get the infinitely 
weak perturbed elliptic point of the map. Because of the elliptic character 
of this fixed point all amplitudes fi have to stay infinitely close to their fixed 
point values which are finite. Thus we have found, that for any initial vector 
fi = (a, b) with a < b there exists exactly one value of k such that the initial 
vector belongs to 5*+. The point fi = (b, a) belongs then to S'_. 

In the case a = 6 it follows that it is always possible to find exactly one 
value of K such that the initial vector // = (a, a) belongs to both S+ and S-. 
Consequently we proved rigorously the existence of one type of NLEs, which 



are known as 'even-parity-mode' 0]. Their characteristic feature is that the 
center of energy density of the corresponding solution is located between two 
lattice sites / and (/ - 1) at (/ - 0.5) [|TU|. 

Finally let us consider the case < fi^2 = fi = ^ < fi-i = b. Solving the 
equation (^ for Z ^ (/ — 1) we get allowed values of k in an interval defined 
by b. Appending the above described procedure for the sequence of lower 
bounds of k we again find that for exactly one value of k such that a < b the 
vectors (a, b) and (6, a) belong to 5*+. By the symmetry of the assumed initial 
configuration they also belong to S^ . Consequently we proved rigorously the 
existence of a second type of NLEs, which are known as 'odd-parity mode' 
[0]. Their characteristic feature is that the center of energy density of the 
solution is located on the lattice site (/ — 1) |TD|. 

Due to the radius of interaction tq = 1 it follows that there are no other 
allowed stationary time-periodic NLE solutions in lattices of type (^. In the 
limit m -^ oo the two NLE solutions become compact and can be calulated 
analytically [|]. 

It would be now very interesting to analyse the map (^ in order to observe 
the previously discussed stable and unstable manifolds, and consequently the 
homoclinic points. However this map (^ is not linearizable around the fixed 
point fi = 0- some matrix elements of the Jacobian diverge. Consequently 
it is impossible to determine the structure of the manifolds near the fixed 
point numerically. In order to proceed we add to the right hand side of (|^) 
the terms 

Cifi + fi-i) + Cifi + fi^,) . (6) 

In this case the modified map becomes linearizable around the fixed point 
fi = 0. We can then visualize the manifolds, and hope that in the limit 
C -^ the structure of the manifolds does not change with respect to their 
intersections. In Fig. 2 we show the stable and unstable manifolds for m = 2, 
C = 2 and k = 10, where we used the variables gi = fi + fi-i instead of 
the original ones. We computed the unstable manifold using standard pro- 
cedures |]T3 and indicated the position of the stable manifold as it can be 



evaluated out of the linearization around the fixed point. Clearly the horse- 
shoe patterns are visible together with the homoclinic points. If one chooses 
one homoclinic point and iterates, then the image is the next-to-next homo- 
clinic point. Thus we indeed observe the two possible even and odd parity 
solutions. If we lower the value of C, then it becomes increasingly harder to 



find the unstable manifold. Still the structure of the unstable manifold can 
be seen partially, which indicates the persistence of the horseshoe patterns in 
phase space regions where the linear terms (|^) are small perturbations. Then 
we can use the results of Section II and conclude that NLE solutions will sur- 
vive under perturbations of the considered system. In particular it follows 
that NLEs are exact solutions for Fermi- Past a- Ulam chains with additional 
linear and cubic spring forces, as often studied in the literature |]ll,[|]. 

4 Conclusion 

In conclusion we have shown, that if a one-dimensional system is known 
where NLE solutions exist and obey the nonresonance condition, then they 
are stable under perturbations of the Hamiltonian of the system or they can 
be perturbed in such a way that they become stable. Consequently we can 
use exact proofs of the existence of NLE solutions, e.g. the one given in 
this work for (Q) or the one in and obtain new NLE solutions for the 
corresponding perturbed system. Remarkably we are not restricted in the 
choice of the perturbation as long as the map for the perturbed system has 
the same phase space dimension. 

Let us discuss the problems of the presented approach which appear when 
we try to consider lattices of dimensionality higher than d = 1. We can still 
make the ansatz of an existing time-periodic NLE solution and will yield 
again a coupled set of algebraic equations for the Fourier coefficients. How- 
ever this set can not be viewed as a discrete map of a certain phase space 
of the Fourier components, ii d > 1. Consequently we do not know how to 
formulate the conditions of the decay of the solution at infinity by imposing 
certain constraints on the set of equations. Still it is well known from nu- 
merical simulations that periodic NLEs exist [15 . 
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FUGURE CAPTIONS 



Fie.l 



Phase space of the map defined by (^ for k = 10 and m = 2. The fixed 
points are: Ajpi = (0,0) (sohd circle), Afp2 = {±x,^x) (x = l/vT6, peri- 
odicity 2), Afp3 = (0, ±y) and {±y, 0) {y = a/5, periodicity 4). 



Fig.2 

Appearance of the homochnic points and horseshoe structure for the per- 
turbed map (H),® with variables as given in the text. 
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